Implementation and Application of Large-Scale Drug Virtual Screening

The molecular docking-based virtual screening technique evaluates the binding abilities between multiple ligand compounds and receptors to screen for the active compounds. In the context of the global spread of the COVID-19 pandemic, large-scale and rapid drug virtual screening is crucial for identifying potential drug molecules from massive datasets of ligand structures. The powerful computing power of supercomputer provides hardware guarantee for drug virtual screening, but the super large-scale drug virtual screening still faces many challenges that affects the effective execution of the calculation. Based on the analysis of the challenges, this paper proposes a centralized task distribution scheme with a central database, and designs a multi-level task distribution framework. The challenges are effectively solved through multi-level intelligent scheduling, multi-level compression processing of massive small molecule files, dynamic load balancing and high error tolerance management technology. An easy-touse"tree”multi-level task distribution system is implemented. A fast, efficient and stable drug virtual screening task distribution, calculation and result analysis function is realized, and the computing efficiency is nearly linear. Then, heterogeneous computing technology is used to complete the drug virtual screening of more than 2 billion compounds, for two different active sites for COVID-19, on the domestic super computing system, which provides a powerful computing guarantee for the super large-scale rapid virtual screening of explosive malignant infectious diseases. © 2023, Journal of Computer Engineering and Applications Beijing Co., Ltd.;Science Press. All rights reserved.

A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years.

Structure-based virtual screening (SBVS), also known as molecular docking, has been increasingly applied to discover small-molecule ligands based on the protein structures in the early stage of drug discovery. In this review, we comprehensively surveyed the prospective applications of molecular docking judged by solid experimental validations in the literature over the past fifteen years. Herein, we systematically analyzed the novelty of the targets and the docking hits, practical protocols of docking screening, and the following experimental validations. Among the 419 case studies we reviewed, most virtual screenings were carried out on widely studied targets, and only 22% were on less-explored new targets. Regarding docking software, GLIDE is the most popular one used in molecular docking, while the DOCK 3 series showed a strong capacity for large-scale virtual screening. Besides, the majority of identified hits are promising in structural novelty and one-quarter of the hits showed better potency than 1 µM, indicating that the primary advantage of SBVS is to discover new chemotypes rather than highly potent compounds. Furthermore, in most studies, only in vitro bioassays were carried out to validate the docking hits, which might limit the further characterization and development of the identified active compounds. Finally, several successful stories of SBVS with extensive experimental validations have been highlighted, which provide unique insights into future SBVS drug discovery campaigns.

Drug Potency Prediction of SARS-CoV-2 Main Protease Inhibitors Based on a Graph Generative Model.

The prediction of a ligand potency to inhibit SARS-CoV-2 main protease (M-pro) would be a highly helpful addition to a virtual screening process. The most potent compounds might then be the focus of further efforts to experimentally validate their potency and improve them. A computational method to predict drug potency, which is based on three main steps, is defined: (1) defining the drug and protein in only one 3D structure; (2) applying graph autoencoder techniques with the aim of generating a latent vector; and (3) using a classical fitting model to the latent vector to predict the potency of the drug. Experiments in a database of 160 drug-M-pro pairs, from which the pIC50 is known, show the ability of our method to predict their drug potency with high accuracy. Moreover, the time spent to compute the pIC50 of the whole database is only some seconds, using a current personal computer. Thus, it can be concluded that a computational tool that predicts, with high reliability, the pIC50 in a cheap and fast way is achieved. This tool, which can be used to prioritize which virtual screening hits, will be further examined in vitro.

Structure-guided pharmacophore based virtual screening, docking, and molecular dynamics to discover repurposed drugs as novel inhibitors against endoribonuclease Nsp15 of SARS-CoV-2.

COVID-19 (Corona Virus Disease of 2019) caused by the novel 'Severe Acute Respiratory Syndrome Coronavirus-2' (SARS-CoV-2) has wreaked havoc on human health and the global economy. As a result, for new medication development, it's critical to investigate possible therapeutic targets against the novel virus. 'Non-structural protein 15' (Nsp15) endonuclease is one of the crucial targets which helps in the replication of virus and virulence in the host immune system. Here, in the current study, we developed the structure-based pharmacophore model based on Nsp15-UMP interactions and virtually screened several databases against the selected model. To validate the screening process, we docked the top hits obtained after secondary filtering (Lipinski's rule of five, ADMET & Topkat) followed by 100 ns molecular dynamics (MD) simulations. Next, to revalidate the MD simulation studies, we have calculated the binding free energy of each complex using the MM-PBSA procedure. The discovered repurposed drugs can aid the rational design of novel inhibitors for Nsp15 of the SARS-CoV-2 enzyme and may be considered for immediate drug development.

Looking for SARS-CoV-2 therapeutics through computational approaches.

BACKGROUND: In the last few years in silico tools, including drug repurposing coupled with structure-based virtual screening, have been extensively employed to look for anti-COVID-19 agents. OBJECTIVE: The present review aims to provide readers with a portrayal of computational approaches that could conduct more quickly and cheaply to novel anti-viral agents. Particular attention is given to docking-based virtual screening. METHOD: The World Health Organization website was consulted to gain the latest information on SARS-CoV-2, its novel variants and their interplay with COVID-19 severity and treatment options. The Protein Data Bank was explored to look for 3D coordinates of SARS-CoV-2 proteins in their free and bound states, in the wild-types and mutated forms. Recent literature related to in silico studies focused on SARS-CoV-2 proteins was searched through PubMed. RESULTS: A large amount of work has been devoted thus far to computationally targeting viral entry and searching for inhibitors of the S-protein/ACE2 receptor complex. Another large area of investigation is linked to in silico identification of molecules able to block viral proteases -including Mpro- thus avoiding maturation of proteins crucial for virus life cycle. Such computational studies have explored the inhibitory potential of the most diverse molecule databases (including plant extracts, dietary compounds, FDA approved drugs). CONCLUSION: More efforts need to be dedicated in the close future to experimentally validate the therapeutic power of in silico identified compounds in order to catch, among the wide ensemble of computational hits, novel therapeutics to prevent and/or treat COVID-19.

In Silico Structure-Based Vaccine Design.

Structure-based vaccine design (SBVD) is an important technique in computational vaccine design that uses structural information on a targeted protein to design novel vaccine candidates. This increasing ability to rapidly model structural information on proteins and antibodies has provided the scientific community with many new vaccine targets and novel opportunities for future vaccine discovery. This chapter provides a comprehensive overview of the status of in silico SBVD and discusses the current challenges and limitations. Key strategies in the field of SBVD are exemplified by a case study on design of COVID-19 vaccines targeting SARS-CoV-2 spike protein.

DHFR Inhibitors Display a Pleiotropic Anti-Viral Activity against SARS-CoV-2: Insights into the Mechanisms of Action.

During the COVID-19 pandemic, drug repurposing represented an effective strategy to obtain quick answers to medical emergencies. Based on previous data on methotrexate (MTX), we evaluated the anti-viral activity of several DHFR inhibitors in two cell lines. We observed that this class of compounds showed a significant influence on the virus-induced cytopathic effect (CPE) partly attributed to the intrinsic anti-metabolic activity of these drugs, but also to a specific anti-viral function. To elucidate the molecular mechanisms, we took advantage of our EXSCALATE platform for in-silico molecular modelling and further validated the influence of these inhibitors on nsp13 and viral entry. Interestingly, pralatrexate and trimetrexate showed superior effects in counteracting the viral infection compared to other DHFR inhibitors. Our results indicate that their higher activity is due to their polypharmacological and pleiotropic profile. These compounds can thus potentially give a clinical advantage in the management of SARS-CoV-2 infection in patients already treated with this class of drugs.

Computational design of medicinal compounds to inhibit RBD-hACE2 interaction in the Omicron variant: unveiling a vulnerable target site.

The COVID-19 pandemic, caused by SARS-CoV-2, has globally affected both human health and economy. Several variants with a high potential for reinfection and the ability to evade immunity were detected shortly after the initial reported case of COVID-19. A total of 30 mutations in the spike protein (S) have been reported in the SARS-CoV-2 (BA.2) variant in India and South Africa, while half of these mutations are in the receptor-binding domain and have spread rapidly throughout the world. Drug repurposing offers potential advantages over the discovery of novel drugs, and one is that it can be delivered quickly without lengthy assessments and time-consuming clinical trials. In this study, computational drug design, such as pharmacophore-based virtual screening and MD simulation has been concentrated, in order to find a novel small molecular inhibitor that prevents hACE2 from binding to the receptor binding domain (RBD). three medicinal compound databases: North-East African, North African, and East African were screened and carried out a multi-step screening approach that identified three compounds, which are thymoquinol 2-O-beta-glucopyranoside (C1), lanneaflavonol (C2), and naringenin-4'-methoxy-7-O-Alpha-L-rhamnoside (C3), with excellent anti-viral properties against the RBD of the omicron variant. Furthermore, PAIN assay interference, computation bioactivity prediction, binding free energy, and dissociation constant were used to validate the top hits, which indicated good antiviral activity. The three compounds that were found may be useful against COVID-19, though more research is required. These findings could aid the development of novel therapeutic drugs against the emerging Omicron variant of SARS-CoV-2.

EXSCALATE: An Extreme-Scale Virtual Screening Platform for Drug Discovery Targeting Polypharmacology to Fight SARS-CoV-2

The social and economic impact of the COVID-19 pandemic demands a reduction of the time required to find a therapeutic cure. In this paper, we describe the EXSCALATE molecular docking platform capable to scale on an entire modern supercomputer for supporting extreme-scale virtual screening campaigns. Such virtual experiments can provide in short time information on which molecules to consider in the next stages of the drug discovery pipeline, and it is a key asset in case of a pandemic. The EXSCALATE platform has been designed to benefit from heterogeneous computation nodes and to reduce scaling issues. In particular, we maximized the accelerators' usage, minimized the communications between nodes, and aggregated the I/O requests to serve them more efficiently. Moreover, we balanced the computation across the nodes by designing an ad-hoc workflow based on the execution time prediction of each molecule. We deployed the platform on two HPC supercomputers, with a combined computational power of 81 PFLOPS, to evaluate the interaction between 70 billion of small molecules and 15 binding-sites of 12 viral proteins of SARS-CoV-2. The experiment lasted 60 hours and it performed more than one trillion ligand-pocket evaluations, setting a new record on the virtual screening scale.

Discovery of dual S-RBD/NRP1-targeting peptides: structure-based virtual screening, synthesis, biological evaluation, and molecular dynamics simulation studies.

Both receptor-binding domain in spike protein (S-RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and human neuropilin-1 (NRP1) are important in the virus entry, and their concomitant inhibition may become a potential strategy against the SARS-CoV-2 infection. Herein, five novel dual S-RBD/NRP1-targeting peptides with nanomolar binding affinities were identified by structure-based virtual screening. Particularly, RN-4 was found to be the most promising peptide targeting S-RBD (Kd = 7.4 ± 0.5 nM) and NRP1-BD (the b1 domain of NRP1) (Kd = 16.1 ± 1.1 nM) proteins. Further evidence in the pseudovirus infection assay showed that RN-4 can significantly inhibit the SARS-CoV-2 pseudovirus entry into 293 T cells (EC50 = 0.39 ± 0.09 µM) without detectable side effects. These results suggest that RN-4, a novel dual S-RBD/NRP1-targeting agent, holds potential as an effective therapeutic to combat the SARS-CoV-2 infection.

Identification of FDA-approved drugs against SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) through computational virtual screening.

The SARS-CoV-2 coronavirus is responsible for the COVID-19 outbreak, which overwhelmed millions of people worldwide; hence, there is an urgency to identify appropriate antiviral drugs. This study focuses on screening compounds that inhibit RNA-dependent RNA-polymerase (RdRp) essential for RNA synthesis required for replication of positive-strand RNA viruses. Computational screening against RdRp using Food and Drug Administration (FDA)-approved drugs identified ten prominent compounds with binding energies of more than - 10.00 kcal/mol, each a potential inhibitor of RdRp. These compounds' binding energy is comparable to known RdRp inhibitors remdesivir (IC50 = 10.09 µM, SI = 4.96) and molnupiravir (EC50 = 0.67 - 2.66 µM) and 0.32-2.03 µM). Remdesivir and molnupiravir have been tested in clinical trial and remain authorized for emergency use in the treatment of COVID-19. In docking simulations, selected compounds are bound to the substrate-binding pocket of RdRp and showed hydrophobic and hydrogen bond interaction. For molecular dynamics simulation, capmatinib, pralsetinib, ponatinib, and tedizolid phosphate were selected from the initial ten candidate compounds. MD simulation indicated that these compounds are stable at 50-ns MD simulation when bound to RdRp protein. The screen hit compounds, remdesivir, molnupiravir, and GS-441524, are bound in the substrate binding pocket with good binding-free energy. As a consequence, capmatinib, pralsetinib, ponatinib, and tedizolid phosphate are potential new inhibitors of RdRp protein with potential of limiting COVID-19 infection by blocking RNA synthesis. Supplementary Information: The online version contains supplementary material available at 10.1007/s11224-022-02072-1.

AI-Driven De Novo Design and Molecular Modeling for Discovery of Small-Molecule Compounds as Potential Drug Candidates Targeting SARS-CoV-2 Main Protease.

Over the past three years, significant progress has been made in the development of novel promising drug candidates against COVID-19. However, SARS-CoV-2 mutations resulting in the emergence of new viral strains that can be resistant to the drugs used currently in the clinic necessitate the development of novel potent and broad therapeutic agents targeting different vulnerable spots of the viral proteins. In this study, two deep learning generative models were developed and used in combination with molecular modeling tools for de novo design of small molecule compounds that can inhibit the catalytic activity of SARS-CoV-2 main protease (Mpro), an enzyme critically important for mediating viral replication and transcription. As a result, the seven best scoring compounds that exhibited low values of binding free energy comparable with those calculated for two potent inhibitors of Mpro, via the same computational protocol, were selected as the most probable inhibitors of the enzyme catalytic site. In light of the data obtained, the identified compounds are assumed to present promising scaffolds for the development of new potent and broad-spectrum drugs inhibiting SARS-CoV-2 Mpro, an attractive therapeutic target for anti-COVID-19 agents.

Identifying potent inhibitory phytocompounds from Lagerstroemia speciosa against SARS-Coronavirus-2: structure-based virtual screening.

The ongoing spillover of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) calls for expedited countermeasure through developing therapeutics from natural reservoirs and/or the use of less time-consuming drug discovery methodologies. This study aims to apply these approaches to identify potential blockers of the virus from the longstanding medicinal herb, Lagerstroemia speciosa, through comprehensive computational-based screening. Nineteen out of 22 L. speciosa phytochemicals were selected on the basis of their pharmacokinetic properties. SARS-CoV-2 Main protease (Mpro), RNA-directed RNA polymerase (RdRp), Envelope viroporin protein (Evp) and receptor-binding domain of Spike glycoprotein (S-RBD), as well as the human receptor Angiotensin-converting enzyme-2 (hACE2) were chosen as targets. The screening was performed by molecular docking, followed by 100-ns molecular dynamic simulations and free energy calculations. 24-Methylene cycloartanol acetate (24MCA) was found as the best inhibitor for both Evp and RdRp, and sitosterol acetate (SA) as the best hit for Mpro, S-RBD and hACE2. Dynamic simulations, binding mode analyses, free energy terms and share of key amino acids in protein-drug interactions confirmed the stable binding of these phytocompounds to the hotspot sites on the target proteins. With their possible multi-targeting capability, the introduced phytoligands might offer promising lead compounds for persistent fight with the rapidly evolving coronavirus. Therefore, experimental verification of their safety and efficacy is recommended.

A Review on Artificial Intelligence Approaches and Rational Approaches in Drug Discovery.

Artificial intelligence (AI) speeds up the drug development process and reduces its time, as well as the cost which is of enormous importance in outbreaks such as COVID-19. It uses a set of machine learning algorithms that collects the available data from resources, categorises, processes and develops novel learning methodologies. Virtual screening is a successful application of AI, which is used in screening huge drug-like databases and filtering to a small number of compounds. The brain's thinking of AI is its neural networking which uses techniques such as Convoluted Neural Network (CNN), Recursive Neural Network (RNN) or Generative Adversial Neural Network (GANN). The application ranges from small molecule drug discovery to the development of vaccines. In the present review article, we discussed various techniques of drug design, structure and ligand-based, pharmacokinetics and toxicity prediction using AI. The rapid phase of discovery is the need of the hour and AI is a targeted approach to achieve this.

An in silico discovery of potential 3CL protease inhibitors of SARS-CoV-2 based upon inactivation of the cysteine 145-Histidine 41 catalytic dyad.

Coronavirus disease 19 (COVID19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Currently, several countries are at risk of the pandemic caused by this virus. In the absence of any vaccine or virus-specific antiviral treatments, the need is to fast track search for potential drug candidates to combat the virus. Though there are known drugs that are being repurposed to fight against the SARS-CoV-2, there is a requirement for the virus-specific drugs at the earliest. One of the main drug targets of SARS-CoV-2 is an essential non-structural protein, 3CL protease, critical for the life cycle of the virus. We have used molecular docking studies to screen a chemically diverse set of small molecules to identify potential drug candidates to target this protein. Of the 22,630 molecules from varied small molecule libraries, based on the binding affinities and physicochemical properties, we finalized 30 molecules to be potential drug candidates. Eight of these molecules bind in a manner allowing for the scope of a nearly orthogonal backside nucleophilic attack on their suitably placed electrophilic carbonyl groups by the thiol group of cysteine residue 145, while remaining inside a 4 Ǻ distance range. It is interesting since carbonyl groups are known to be attacked in a similar fashion by external nucleophiles and can be relevant when considering these molecules as potential mechanism-based irreversible inhibitors of the 3CLPro. Further, ADMET analysis and Molecular dynamics simulations and available bioactive assays led to the identification of three molecules with high potential to be explored as drug candidates/lead molecules to target 3CLPro of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Poster: An Extreme-Scale Virtual Screening Platform for Drug Discovery

Virtual screening is one of the early stages that aims to select a set of promising ligands from a vast chemical library. Molecular Docking is a crucial task in the process of drug discovery and it consists of the estimation of the position of a molecule inside the docking site. In the contest of urgent computing, we designed from scratch the EXSCALATE molecular docking platform to benefit from heterogeneous computation nodes and to avoid scaling issues. This poster presents the achievements and ongoing development of the EXSCALATE platform, together with an example of usage in the context of the COVID-19 pandemic.

Virtual Screening for Sars-Cov-2 Entry Inhibitors by Dual Targeting of Tmprss2 and Ctsl

The COVID-19 pandemic remains to be a global public health crisis due to the emergence of new variants of concern and the scarcity of drug treatments. The cell entry of SARS-CoV-2 requires activation of its spike protein by host proteases TMPRSS2 and CTSL, which triggers membrane fusion and facilitates the endocytic uptake mechanism, respectively. This study employed a structure-based virtual screening technique to identify drugs and natural products that simultaneously target TMPRSS2 and CTSL. Two pharmacophore models were generated from the binding sites of the proteins in complex with their co-crystallized ligands. Both structure-based pharmacophores were used to screen a ligand library composed of 41,775 compounds (10,849 drugs from the ChEMBL database and 30,926 natural products from the NPASS database). A total of 115 compounds (54 drugs and 61 natural products) that fit both TMPRSS2 and CTSL pharmacophore models were identified. The common hits were docked into both proteases to obtain a short list of compounds. Molecular docking filtered 17 compounds (5 drugs and 12 natural products) that have higher binding energy values than the co-crystallized ligands and known inhibitors of both proteins. The top hits were then subjected to ADMET, drug-likeness, and synthetic accessibility filters. Based on docking scores, pharmacokinetics, and drug-likeness, Silibinin was the most promising repurposed drug candidate as a treatment for SARS-CoV-2 infection via dual inhibition of TMPRSS2 and CTSL. Among the natural products, barettin was the best candidate for further development as a novel dual TMPRSS2 and CTSL inhibitor.

Virtual Screening of Natural Compounds Against Six Protein Receptors Coded by The SARS-CoV-2 Genome

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the virus that causes Coronavirus 2019 (COVID-19). To date, there has been no proven effective drug for the treatment or prevention of COVID-19. A study on developing inhibitors for this virus is carried out using molecular docking simulation methods. 3CL-Pro, PL-Pro, Helicase, N, E, and M protein were used as protein targets. Autodock Vina, Autodock 4.2, and PSOVina were used in this study. This study aims to obtain a model of ligands interactions of active natural compounds against the receptor protein encoded by the SARS-CoV-2 genome and their free binding energy to propose active compounds from natural products that have potential as a drug for COVID-19. Corilagin (-14,42 kcal/mol), Scutellarein 7-rutinoside (-13,2 kcal/mol), Genistein 7-O-glucuronide (-10,52 kcal/mol), Biflavonoid-flavone base + 3O (-11,88 and-9,61 kcal/mol), and Enoxolone (-6,96 kcal/mol) has the best free energy value at each protein target indicating that the compound has the potential as a viral protein inhibitor for further investigation. This research is limited to computer simulations, where the results obtained are still a prediction. © 2023, Universitas Jenderal Soedirman. All rights reserved.

Identification of Synthetic 2-Mercaptobenzimidazole Derivatives as Inhibitors of Spike Protein of SARS-CoV-2 by Virtual Screening

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spread worldwide and caused the COVID-19 pandemic. Despite countless efforts in searching for repositioned drugs to treat this disease, the results are still modest. Thus, searching for new compounds as promising drugs to treat this disease is crucial. 2-Mercaptobenzimidazole, a scaffold found in many biologically relevant compounds, has been extensively studied due to its range of biological activities. Using in silico tools, this study aimed to identify 2-mercaptobenzimidazole derivatives as potential drugs to inhibit SARS-CoV-2 infection by blocking spike protein – human angiotensin-converting (hACE2) enzyme interaction. 61 compounds were screened to evaluate their absorption, distribution, metabolism, and excretion (ADME) properties. The compounds that did not violate any Lipinski or Veber rules were subjected to molecular docking interactions to verify their ability to inhibit spike glycoprotein from binding to the hACE2 enzyme. The docking scores for these compounds were superior to the antiviral drugs Remdesivir and Umifenovir, proven to be potential drugs against COVID-19. Furthermore, 2-mercaptobenzimidazole derivatives have been shown to interfere in amino acids involved in key contact sites between SARS-CoV-2-CTD and hACE2 subdomain I. Finally, some toxicological properties were predicted, and the compounds exhibited few (or none) alerts for toxic endpoints as well as low predicted acute oral toxicity (LD50). In this way, the results presented in this work should contribute to discovering new drugs against SARS-CoV-2. © 2022 by the authors.

Predictive analysis, diagnosis of COVID-19 through computational screening and validation with spectro photometrical approach

Objective: To develop Favipiravir, based predictive models of coronavirus disease 2019 (COVID-19) from small molecule databases such as PubChem, Drug Bank, Zinc Database, and literature. Method(s): High Throughput Virtual Screening (HTVS) using different computational screening methods is used to identify the target and lead molecules. CoMFA (Comparative Molecular Field Analysis) is a 3D-QSAR procedure depending on information from known dynamic atoms and eventually permits one to plan and anticipate exercises of particles. These two analysis is used to train predictive models. Result(s): The predictive model achieved the highest accuracy score with a relatively small dataset size can be a subject of overfitting. Datasets with over 500 samples demonstrate an accuracy of about 85-95%, that can be considered as very good. Conclusion(s): From the result it is observed that Increasing level of potassium, sodium and nitrogen will lead to burst lipid bilayer membrane of virus which cause RNA replication rapidly. However, low level of sodium, potassium and nitrogen will help in the DNA polymerase inhibition and replication can be stopped. The best developed QSAR model in terms of the druggability and activity relation has been selected over the parent Favipiravir molecule for designing COVID-19 drugs may lead towards pharmaceutical development in future.Copyright © 2023, The Author(s), under exclusive licence to Korean Society of Environmental Risk Assessment and Health Science.